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PREFACE 


The  theory  of  the  propagation  of  a  finite  amplitude  wave  in 
isotropic  solids  and  crystals  of  cubic  symmetry  was  begun  a  number 
of  years  ago  [M.  A.  Breazeale  and  Joseph  Ford,  J.  Appl.  Phys.  36,  3486- 
3490  (1965)],  and  subsequently  was  developed  In  terms  of  the  elegant 
Lagranglan  formalism  [Albert  C.  Holt  and  Joseph  Ford,  J.  Appl.  Phys. 

40,  142-148  (1969)].  This  theory  has  served  as  the  basis  of  measure¬ 
ment  of  the  third-order  elastic  (TOE)  constants  of  many  isotropic 
solids  (e.g.,  John  H.  Cantrell,  Jr.,  Technical  Report  No.  14)  and 
crystals  of  cubic  symmetry  (e.g.,  Jacob  Philip.  Technical  Report  No.  19). 
For  isotropic  solids  three  TOE  constants  are  necessary,  and  for  cubic 
crystals  one  finds  six.  However,  many  materials  of  technological 
importance  have  lower  symmetry  than  cubic.  Therefore,  it  now  is 
desirable  to  determine  which  TOE  constants  enter  into  harmonic  generation 
in  waves  of  finite  amplitude  propagating  in  crystals  of  lower  symmetry. 
This  Technical  Report  Is  the  first  effort  in  this  direction. 

— The  Technical  Report  gives  a  complete  analysis  of  the  theory  of 

.  r  r e~a-\r~ 

propagation  of  a  finite  amplitude  ultrasonic  wave  in  ancrystal  of 
hexagonal  symmetry  along  the  symmetry  axis  and  along  any  direction  in 
the  basal  plane.  From  the  theory  we  conclude  that  at  least  four  of  the 

h',  r  >  e^hySt'  C 

ten  TOfr  constants  can  be  evaluated  from  harmonic  generation  measurements 

in  the  appropriate  directions  In  hexagonal  crystals.  <^~ - — — 

Completion  of  the  Technical  Report  was  made  possible  by  the 
efforts  of  Gonghuan  Du  whose  meticulous  mathematical  efforts  are 
gratefully  acknowledged. 


CHAPTER  I 


INTRODUCTION 

Although  the  nonlinear  acoustical  properties  of  fluids  has  been 
studied  for  more  than  a  century,  a  systematic  study  of  the  nonlinear 
acoustical  properties  of  solids  actually  began  about  three  decades 
ago.  Work  in  this  area  gained  momentum  rapidly  because  a  large  number 
of  the  basic  properties  of  solids  can  be  explained  only  on  the  basis 
of  a  nonlinear  theory.  Among  these  properties  are  thermal  expansion, 
the  Interaction  of  lattice  waves.  Inequality  of  adiabatic  and  isothermal 
elastic  constants,  and  the  dependence  of  the  elastic  constants  on 
pressure  and  temperature. 

As  Is  well  known,  the  description  of  the  nonlinear  properties 

of  solids  Is  made  In  terms  of  the  higher-order  elastic  constants.  These 

constants,  coefficients  of  the  higher-order  terms  in  the  relation  between 

Internal  energy  and  latt^e  strains  on  a  solid,  are  some  of  the  most 

important  parameters  In  the  description  of  the  nonlinear  properties  of 

solids.  Their  measurement,  therefore,  has  attracted  considerable 

attention.  The  ordinary  elastic  constants,  of  course,  are  determined 

directly  from  the  sound  velocity.  They  are  the  second-order  elastic 

constants.  The  third-order  elastic  constants  are  measured  by  techniques 

outlined  In  a  previous  Technical  Report.^  As  was  pointed  out  in  that 

Report,  most  of  the  measurements  reported  are  on  solids  of  cubic 

symmetry  with  six  third-order  elastic  (TOE)  constants..  As  the  symmetry 

? 

Is  lowered,  the  number  of  elastic  constants  Increases.  For  example, 

1 


2 


crystals  with  hexagonal  symmetry  have  10  TOE  constants  and  those  with 
trigonal  symmetry  have  14  (or  20)  TOE  constants. 

Only  very  few  measurements  have  been  reported  on  the  TOE  constants 
of  hexagonal  and  trigonal  crystals,  in  part  because  of  a  lack  of  a 
general  theory.  The  linear  theory  of  elastic  wave  propagation  along 

O 

different  directions  In  these  types  of  media  has  been  given  by  Musgrave, 

4 

Farnell,  and  others.  A  complete  description  of  the  linear  theory  has 

c 

been  given  by  Musgrave.  The  purpose  of  the  present  study  is  to  expand 
on  the  linear  theory  to  Include  nonlinear  terms  and  to  formulate  the 
theory  in  such  a  way  that  the  nonlinear  distortion  of  finite  amplitude 
ultrasonic  waves  in  hexagonal  crystals  can  be  described. 

In  an  anisotropic  medium  there  are  only  certain  directions  along 
which  elastic  waves  can  propagate  as  the  pure  longitudinal  modes  which 
are  of  primary  interest  In  the  application  of  the  harmonic  generation 
technique.  Associated  with  any  propagation  direction  there  are  three 
independent  waves,  the  displacements  of  which  form  a  mutually  orthogonal 
set.  In  general,  none  of  the  three  displacement  vectors  coincides 
with  the  vector  which  Is  normal  to  the  wavefront;  i.e.,  in  general,  the 
waves  are  neither  pure  longitudinal  nor  pure  transverse.  The  specific 
directions  in  which  pure  mode  longitudinal  waves  propagate  have  been 
determined  by  Borgnls®  for  crystals  with  trigonal,  hexagonal,  tetragonal 
and  cubic  symmetry  and  general  conditions  have  been  established  under 
which  pure  longitudinal  waves  exist.  In  order  to  determine  the  third- 
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order  elastic  constants  from  wave  speed  measurements  In  stressed 
crystals.  It  is  necessary  to  know  these  pure  mode  directions. 

The  theory  of  harmonic  generation  of  longitudinal  waves,  which 
is  discussed  in  detail  in  later  chapters,  also  depends  upon  availability 
of  pure  longitudinal  mode  directions.  Using  the  method  due  to 
Borgnis,®  Brugger^  has  determined  the  pure  mode  directions  for  all 
crystal  point  groups  belonging  to  orthorhombic,  tetragonal,  cubic, 
trigonal  (rhombohedral )  and  hexagonal  systems.  Also  Brugger^  has 
solved  the  eigenvalue  problem  for  each  of  these  directions  and 
polarization  vectors  and  the  wave  speeds  have  been  tabulated.  A  similar 

p 

description  for  transverse  waves  has  been  given  by  Chang. 

The  theory  of  finite  deformations  and  its  application  to 

crystalline  solids  Including  nonlinear  terms  has  been  given  by 
g 

Murnaghan.  The  theory  has  been  applied  to  cubic  crystals  subjected 

to  finite  hydrostatic  compression  by  Birch10  and  later  by  Seeger  and 
11  12 

Buck.  Einspruch  and  Manning  extended  the  theory  to  hexagonal, 

13  14 

tetragonal  and  orthorhombic  symmetry.  Thurston  and  Brugger  *  solved 
the  equations  of  small  amplitude  waves  in  a  homogeneously  deformed 
crystal  and  expressed  their  results  in  terms  of  the  "natural  velocity" 
and  its  stress  derivative.  They  defined  the  "natural  velocity"  as  the 
length  of  the  specimen  In  the  unstressed  state  divided  by  the  wave 
travel  time  in  the  stressed  state.  They  presented  results  for  isotropic 
and  cubic  solids.  The  necessary  relations  between  third-order  elastic 
coefficients  and  stress  derivatives  of  sound  speeds  along  different 
directions  in  crystals  of  orthorhombic,  tetragonal,  cubic,  trigonal  and 


4 


1 5 

hexagonal  classes  have  been  derived  by  Brugger  using  his  thermodynamic 

definition  of  third-order  elastic  constants.16 

The  TOE  constants  of  a  large  number  of  crystals  belonging  to  cubic 

symnetry  have  been  measured.  Most  of  these  measurements  involve 

measurement  of  the  stress  derivatives  of  ultrasonic  wave  velocities  in 

solids.  The  first  measurement  of  this  kind  was  by  Lazarus17  on  KC1, 

18 

NaCI,  CuZn,  Cu  and  A1.  Later  Hughes  and  Kelly  made  measurements  on 

polystyrene,  iron  and  pyrex  glass  and  reported  the  three  Murnaghan  TOE 

constants  of  these  materials.  The  first  measurement  of  the  complete  set 

19 

of  six  TOE  constants  of  a  cubic  solid  was  by  Bateman  et  al.  on  ger¬ 
manium.  Later,  measurements  have  been  reported  on  many  cubic  crystals 

by  several  workers.  The  values  reported  up  to  1979  have  been  tabulated 
2 

by  Hearmon. 

The  second  most  widely  used  technique  to  measure  the  TOE  constants 

of  solids  is  the  ultrasonic  harmonic  generation  technique  developed  by 

20 

Breazeale  and  Ford.  Measurements  have  been  reported  on  a  number  of 
21 

solids.  The  theory  of  second  harmonic  generation  in  cubic  crystals 
and  the  experimental  technique  has  been  described  in  detail  in  a  pre¬ 
vious  report.1  Also  for  a  review  of  the  work  done  on  the  TOE  constants 

2  22 

of  cubic  crystals,  refer  to  Hearmon  and  Green. 

Even  though  the  TOE  constants  of  a  number  of  isotropic  and  cubic 
crystals  have  been  determined,  not  many  measurements  have  been  reported 
on  crystals  of  lower  symmetry.  In  all  the  work  reported  on  such 
crystals  so  far,  the  hydrostatic  and  uniaxial  stress  derivatives  of 
ultrasonic  wave  velocities  have  been  measured  and  the  TOE  constants 
have  been  determined  from  those  data.  Since  the  number  of  TOE  constants 


5 


for  crystals  of  lower  symmetry  is  higher,  a  correspondingly  higher 
number  of  careful  measurements  are  to  be  made  to  isolate  all  the  TOE 
constants.  In  the  following  paragraphs  we  outline  the  measurements 
to  be  made  on  hexagonal  crystals  with  10  TOE  constants.  The  values 

measured  until  197?  on  these  solids  also  have  been  tabulated  by 

2 

Hearmon. 

The  complete  set  of  ten  TOE  constants  of  the  hexagonal  crystal 
zinc  has  been  determined  from  measurements  of  the  hydrostatic  and  uni¬ 
axial  stress  dependence  of  ultrasonic  wave  velocity  by  Swartz  and 
23 

Elbaum.  They  used  data  from  twelve  measurements  for  various  propagation, 
polarization  and  stress  application  directions  to  isolate  all  the  10  TOE 
constants.  An  interference  technique  has  been  used  by  them  to  measure 
the  wave  velocity  difference  in  stressed  and  unstressed  samples.  They 
have  reported  the  temperature  and  pressure  derivatives  of  the  SOE  con¬ 
stants  in  their  work.  The  Griinelsen  parameter  y  of  zinc  has  been 
evaluated  by  the  authors  and  the  volume  thermal  expansion  coefficient 
calculated  from  that  is  in  agreement  with  experimental  value  of  the 
volume  thermal  expansion  coefficient. 

The  complete  set  of  the  ten  TOE  constants  of  the  hexagonal  metal 

24 

magnesium  has  been  determined  by  Naimon,  who  measured  the  hydrostatic 
pressure  and  uniaxial  compression  derivatives  of  natural  sound  velocities 
using  ultrasonic  pulse-superposition  technique.  With  the  sample 
oriented  with  faces  perpendicular  to  the  a,  b,  c  axes,  four  hydrostatic 
and  ten  uniaxial  experiments  for  different  propagation  and  polarization 
directions  were  used  to  determine  the  ten  TOE  constants.  Even  though 

25 

the  hydrostatic  pressure  derivatives  of  the  SOE  constants  of  beryllium. 


zirconium,  europium,  rhenium,  magnesium,  cadmium,  , 

37  37  37 

gadolenium,  dysprosium,  and  titanium  have  been  reported,  that 

information  is  not  sufficient  to  determine  the  complete  set  of  TOE 

constants.  The  TOE  constants  of  a  number  of  hexagonal  solids  have 

been  determined  theoretically  by  a  number  of  workers  using  the  method 

of  homogeneous  deformations  with  different  model  potentials.  For  a 

32 

review  of  the  theoretical  work,  refer  to  Ramji  Rao  and  Ramanand. 

In  the  following  report  we  develop  the  nonlinear  theory  required 
to  consistently  describe  the  propagation  of  an  initially  sinusoidal 
ultrasonic  wave  along  pure  mode  directions  in  a  solid  of  hexagonal 
symmetry.  As  pointed  out  by  Brugger,^  pure  mode  directions  in  a 
hexagonal  crystal  are  (1)  the  symmetry  axis,  (2)  all  directions  in  the 
basal  plane,  and  (3)  directions  along  a  cone  centered  on  the  symmetry 
axis,  the  apex  angle  of  the  cone  depending  on  the  SOE  constants  of  the 
hexagonal  crystal  under  consideration.  In  the  following  report  we 
specialize  the  equations  to  the  pure  mode  directions  (1)  and  (2)  and 
sK"  which  TOE  constants  can  be  measured  by  propagating  finite  amplitude 
waves  along  these  directions.  We  end  the  report  with  some  comments 
about  the  way  in  which  the  theory  can  be  handled  for  directions  (3) 
along  a  cone  centered  on  the  symmetry  axis. 


CHAPTER  II 


GENERAL  THEORY  OF  NONLINEAR  WAVE  PROPAGATION  IN  SOLIDS 

Consider  a  point  P  in  the  medium  with  coordinates  a^  (a,b,c)  in 
the  unstrained  state.  Let  the  point  P  move  to  P*  with  coordinates 
xi  (x,y,z)  in  the  deformed  state.  The  components  of  the  displacement 
can  then  be  written  as 

U  =  x  -  a 

V  =  y  -  b  (1) 

W  =  z  -  c 

In  the  theory  of  finite  deformations  of  an  elastic  solid  due  to 
the  large  deformations  involved,  the  initial  coordinates  of  a  particle 
in  the  undeformed  state  are  not  Interchangeable  with  the  final 
coordinates  in  the  deformed  state.  The  expression  for  the  strain 
energy  in  terms  of  the  strains  change  correspondingly.  In  the 
Lagrangian  formulation,  the  strain  is  described  in  the  initial  or 
undeformed  state  and  the  initial  coordinates  of  the  material  particle 
a^  are  taken  as  independent  variables.  The  Lagrangian  formulation  is 
used  in  the  theory  described  in  this  technical  report. 

The  Lagrangian  strain  parameters  which  are  components  of  the 

g 

finite  strain  tensor  are  given  by 

n  *  \  (J*J  -  6)  (2) 

where  J  is  the  Jacobian  matrix  given  by 
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is  the  unit  matrix.  The  stress-strain  relation  is  given  by" 


(4) 


(5) 


(6) 


^p/p0)  0  (U)J* 


(7) 


where  o  is  the  Cauchy  stress  tensor  and  pg  and  p  are  the  densities  in  the 
undeformed  and  deformed  state,  respectively.  <(>  is  the  strain  energy  per 
unit  undeformed  volume  and 
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The  properties  of  the  crystalline  medium  enter  into  the  theory  through 
the  strain  energy  term  $.  Many  deformable  media  with  crystalline 
structure  are  elastically  insensitive  to  certain  prerotations  depending 
upon  the  symmetry  of  the  medium.  For  those  rotations  to  which  the 
medium  is  insensitive  we  have  the  relation 

<t>(R*nR)  ■  4>(n)  (9) 

where  R  is  the  rotation  matrix  and  R*  is  the  transpose  of  R.  can  be 
expanded  and  written  as  a  sum  of  terms  of  different  degrees  in  the 
elements  of  n  as  follows: 

<t>  s  4>q  +  4>i  +  $2  +  $3  +  •••  (10) 

In  that  case  the  series  relation 

4>j(R*nR)  =  4>j(n)  (11) 

holds  good,  $  being  homogeneous  functions  of  degree  j  in  the  elements 
of  n,  j  *  0,1 ,2,3  ...  In  terms  of  the  elastic  moduli,  Eq.  (10)  can  be 
written  as: 


10 


♦  *  *0  +  kl  Cij  nij  +  k2  C1  jk£  nij  \ i 

+  k3  Ci jkumn  n1j  nk£  nmn  +  ”*  (12* 

where  the  C's  are  the  elastic  constants  and  k„  is  a  constant  factor 
depending  on  the  definition  of  elastic  constants.  If  the  initial 
energy  and  deformation  of  the  body  are  neglected,  the  first  two  terms 
in  (12)  are  zero,  and 

♦  =  k2  Ci j k£  nij  nk£  +  k3  Cijk£mn  nij  nkz  "mn  +  ***  (13* 

As  is  well  known,  C^jka  are  the  second-order  and  CijkJJT)n  are  the  third- 

order  elastic  constants.  In  accordance  with  Brugger's  thermodynamic 
definition  of  elastic  constants,16  the  factor  kn  Is  pjy  .  So,  based  on 
that  definition  we  have  k2  *  ^  and  k3  *  Using  single  subscript 
notation  instead  of  double  subscript  notation  for  elastic  constants  and 
using  Brugger's  definition  of  elastic  constants,  the  general  expression 
for  $  can  be  written  as 

♦  ■  1  cx„  "u  "kt +  r  ck„v "«  "n  v,  ♦  h1«her  or,ler  te™s  <'4> 

where  x  ■*  ij,  y  -*•  kt  and  v  -*•  mn.  The  second-order  stiffness  CA^  form  a 

fourth  rank  tensor  containing  81  components,  of  which  21  are  Independent 
for  the  most  unsymmetrlc  trl clinic  and  the  third-order  stiffness  form 
a  sixth  rank  tensor  with  729  components  of  which  56  are  independent  for 
a  trlcllnlc  crystal.  The  number  of  elastic  constants  decreases  con¬ 
siderably  for  crystals  of  higher  symmetry. 
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Equation  (14)  for  the  strain  energy  takes  appropriate  forms  for 

crystals  of  different  synmetry.  The  number  of  elastic  constants  in 

different  crystal  classes  have  been  worked  out  by  various  authors  and 

2 

have  been  presented  in  a  tabular  form  by  Hearmon  which  in  conjunction 
with  Eq.  (14)  conveniently  can  be  used  to  obtain  an  expression  for  the 
strain  energy  of  any  class  of  crystal. 

The  principle  of  conservation  of  mass  along  with  the  definition 
of  J  leads  to  the  relation 


|J|  ■  p0/p  -  (1  +  Ua  +  Vb  +  Wc)  .  (15) 

The  equations  of  motion  for  an  elastic  medium  are  restatements  of 
Newton's  second  law.  For  convenience  one  can  introduce  the  stress 
tensor  T,  which  is  not  symmetric,  as 


T  - 


(16) 


with  nine  components.  This  expression  allows  one  to  write  the  equations 
of  motion  as 


=  p0  Ui 


(17) 


in  the  Lagrangian  coordinate  frame.*  These  equations  of  motion  take 


♦In  the  Eulerian  coordinate  system  the  corresponding  equations  of 
motion  are  3o.,. 

157spUj 

J  ‘  '  ‘  *  the 


'1 


where  is  the  Cauchy  stress  tensor  which  is  symmetric  and  x.  are 
coordinates  of  the  particle  a  in  the  deformed  state. 
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I 


particular  forms  along  the  axes.  Along  the  a,  b  and  c  axes  they  are 


pQ  U  =  3T^/3a  +  3T-j2/3b  +  ST-jj/sc 


V  =  3T01/3a  +  3T99/3b  +  31W3C 


pQ  W  =  3T31/3a  +  3T32/3b  +  3T33/3c 


From  Eqs.  (4),  (8)  and  (16),  the  stress  matrix  T  can  be  written  in 
component  form  as 

T11  T12  T13  J11  J12  J13  9^  9^J  9^ 

T21  T22  T23  *  J21  J22  °23  9^  9^  ^  (19> 

t  t  T  ill  34  34  34 

T31  *32  1 33 J  [J31  J32  J33j[3*31  3*32  3n33. 

From  this  we  can  write  down  the  components  as  follows: 


T 

s 

J,, 

34 

+  J,„ 

34 

+  J,_ 

34 

'll 

dll 

3*n 

d12 

3*21 

13 

3n31 

8 

J,. 

34 

+  J.. 

34 

+  J,. 

34 

12 

dll 

99,2 

u12 

3n22 

u13 

3*32 

T _ 

* 

J,, 

99 

+  J, . 

34 

13 

dll 

3n13 

d12 

5*23 

13 

3n33 

T 

s 

34 

+  J„,» 

34 

+  J„- 

34 

'21 

J21 

3nll 

d22 

3*21 

23 

3n31 

T _ 

e 

3<f> 

+  J— 

34 

+  J„- 

'22 

J21 

99,2 

°22 

3*22 

23 

3*32 

s 

99 

+  J.. 

34 

+  Jo- 

34 

'23 

d21 

3n13 

d22 

3*23 

23 

3n33 

T 

8 

,1 

3(j> 

34 

+  J.. 

34 

'31 

J31 

3nll 

J32 

9n2, 

d33 

3n31 

T„_ 

8 

34 

34 

+  J.. 

'32 

U31 

3n12 

J32 

3*22 

33 

3*32 

T  , 

8 

_il_ 

+  J_- 

3$ 

+  Jo- 

'33 

d31 

3t>13 

g32 

3*23 

33 

3n33 

t 


Note  that  the  T^j  tensor  defined  here  Is  not  symmetric. 

We  will  consider  the  case  of  plane  finite  amplitude  waves 
propagating  along  the  axes  of  the  medium  under  consideration.  For  plane 
waves  propagating  along  the  a-axls  the  displacement  components  become 

U  *  U (a , t) 

V  *  V(a,t)  (21) 

W  *  W(a,t) 

For  this  special  case,  the  equations  of  motion  given  by  Eq.  (18)  become 

Pq  U  ®  3T-j-j/3a 

p0  V  *  3T21/3a  (22) 

p0  =  sT31/aa 

Similarly  for  plane  waves  propagating  along  the  b-axis,  we  have 
P0  ii  *  3T12/3b 

Pg  V  *  3T22/3b  (23) 

p0  W  *  3T32/3b 

and  for  plane  waves  propagating  along  the  c-axis,  we  have 

Pg  U  =  3T^3/3C 
p0  V  ■  3T23/3c 


Pn  W  =  3T„/3C 


These  equations  are  written  out  In  more  explicit  terms  and  solved 
for  crystals  belonging  to  hexagonal  symmetry  in  the  following  chapters. 
The  elastic  strain  energy  appropriate  to  the  crystal  symmetry  is 
obtained.  Only  terms  up  to  third  order  are  retained  in  the  strain 
energy  expression  as  we  are  interested  only  in  the  first  order  of 
nonlinearity  which  is  enough  to  account  for  most  of  the  nonlinear 
properties.  Moreover,  inclusion  of  fourth-  and  higher-order  terms  make 
the  algebra  unwieldy. 


CHAPTER  III 


SECOND  HARMONIC  GENERATION  OF  ULTRASOUND  IN  CRYSTALS 
OF  HEXAGONAL  SYMMETRY 

There  are  two  classes  of  crystals  with  hexagonal  symmetry,  one 
with  the  Hermann-Mauguin  symbol  6,  £,  6/m  with  5  second-order  and  12 
third-order  elastic  constants  and  the  other  with  the  symbol  622,  6mm, 
Fm2,  6/mmm  with  5  second-order  and  10  third-order  elastic  constants. 
Since  most  of  the  crystals  with  hexagonal  symmetry  belong  to  the 
second  class,  we  will  confine  our  attention  to  it.  The  figure  of  the 
hexagonal  close  packed  structure  is  shown  in  Fig.  1.  A  drawing  of 
the  coordinate  frame  and  the  basal  plane  are  also  given  In  the 
figure.  The  five  SOE  constants  are  C^,  C12»  C13,  C33  and  C44  and  the 
ten  TOE  constants  are  C^^,  C^,  C113,  C123,  C133,  C144,  C155,  C222, 

C333  and  C344’2 

As  has  been  shown  by  Borgnis®  and  by  Brugger,^  any  axis  in  the 
xy  plane  and  the  z-axls  are  directions  in  which  pure  mode  longitudinal 
wave  propagation  Is  possible  (In  the  linear  approximation).  To  be 
consistent  with  the  general  theory  in  Chapter  II,  the  x-di recti  on  is 
represented  as  a-axis  and  the  60°  rotated  direction  from  the  x-direction 
is  designated  as  a '-axis.  The  z-di recti  on  and  c-axis  coincide. 

Using  Brugger's  definition  of  third  order  elastic  constants, 
expression  for  the  strain  energy  density  for  hexagonal  crystals  Is 
given  by: 
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Hexagonal 

Structure 


A 


a  =  b  *  c 
a  =  £=90° 


17 


* s  7cn(nii2  +  T1222  +  nU  +  n212) 


+  C12(nnn22  -  |  n122  -  \  n212) 

+  C13(nlln33  +  n22n33^  +  7  C33n332 
+  C44(n132  +  n312  +  n232  +  n322) 


+  6  Cm(t1ll  +  n22^3 


1  2  1 
+  7  C113T,33^nll  +  n22^  +  7  C133T1332(n„  +  %) 


11  T  22 

2  2  2  2 
+  ci44(r’n  +  ^22^23  +  n32  +  n13  +  n31  ^ 

\2  ,  «*/_  2  .  2 


+  7  C166^nll  +  n22^^11  "  n22^  +  Z^n12  +  n21  ^ 

+  lC266^nll  '  n22^3  "  6^n122  +  n212^nll  "  n22^ 


,  1  r  _  3  _  ,  2  ,  2.  2. 

5  C333n33  +  C344T133  n23  n32  n13  n31  ' 

+  7  C366n33^nll  ”  n22^  +  2^n122  +  n212^ 

2  2  2  2 

+  2C456^n22  "  nll^n23  +  n32  "  n13  “  n31  ^  +  4^n12n23n31 
+  n21n32T)13^  • 


(25) 


In  this  expression  C^gg,  C2gg,  C3gg  and  C4gg  are  combinations  of  TOE 
constants  given  by 

C1 66  *  i("2Clll  '  C112  *  3C222} 

C266  *  T<2C111  “  C112  “  C222} 

C366  *  7^113  ’  C123} 

C456  *  7^"C144  +  C155J 


(26) 


,  '  .  X  m  ii  , 
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Elnspruch  and  Manning  derived  the  above  expression  using  Birch's 
definition  of  TOE  constants.  That  expression  can  easily  be  converted 
into  Eq.  (25)  which  was  Brugger's  definition.  The  relation  between 
Brlch  and  Brugger  TOE  constants  can  be  obtained  from  the  relation 


c(Birch)  s  M  c(Brugger)/g 

Apv  Ayv 


(27) 


where  M  is  the  possible  number  of  ways  in  which  can  be  expressed 
In  tensor  notation. 

The  strain  derivatives  of  the  elastic  energy  can  be  obtained  by 
differentiating  Eq.  (25)  with  respect  to  the  nine  components  of  the 
strain  parameter.  They  are  given  by 

3n^  *  Cllnll  +  C12n22  +  C13n33  +?  Clllnll2  +  Clllnlln22  +  7  Cllln222 
1  2 

+  1  C133n33  +  C113nlln33  +  C113n22n33 

2  2  2  2  3  2 

+  C144n23  +  C144n32  +  C144n13  +  C144n31  +  ?  C166r’ll 

+  lC166n222  "  C1 66nl  1  n22  +  C166n122  +  C166n212 

.  1  r  2r  .  1  r  2  r  12 

2  U266nl  1  "  t'266nlln22  1  u266n22  "  l'266T,12 

2  2  2 
‘  C266n21  +  C366nllT|33  "  C366n22n33  "  2C456n23  “  2C456n32 


'366  1 r  33  v366"22"33 


+  2C456n132  +  2C456n312 


*  C,in0i  -  C10noi  +  2C1K(;Tinn9i  +  2C1i;£n99n9i  -  ^^266T1lln21 


(28a) 


3n21  '*'1 1  n21  12  21  tv'166"ir21  ‘-‘'166"22"21 

+  2C266n22n21  +  2C366n33n21  +  SC456r|32nl  3  * 


(28b) 


m 


-■  T1  ...  M.  ■■  •-  - 


44n31  *h44nlln31  ^144n22n31  £L344n33n31  “  HV<456n22 '31 


4C456nlln31  +  8C456n12n23  * 


(28c) 


ln12  *  C12n12  +  2C166nnn12  +  2C166n22n12  ’  2C266nlln12 


2C266n22n12  +  2C366n33nl 2  +  8C456n23n31 


(28d) 


1  2 

in22  +  ci2nn  +  ci3n33 +  2  cmnn  +  Clllnlln22 


1  2  1  2 

7  Cmr,22  *7  C133n33  +  C113n11n33  +  Cn3n22n33 

2  2  2  2 
C144n23  +  C144n32  +  C144n13  +  C144n31 

1  2  3  2  2 

2  C1 66nl 1  +  2  C166n22  "  C166nlln22  +  C166n12 

2  12  12 
C166n21  "  7.  C266nll  +  C266nlln22  "  7  C266n22 

2  2  2 
C266n12  +  C266n21  +  C366n22n33  *  C366nlln33  +  2 C456n23 


2  2  2 
2  C456n32  "  2  C456nl 3  "  2 C456n31 


(28e) 


'44n32  +  2C144nl 1 n32  +  2C1 44n22n32  +  2C344n33n32  +  4C456n22n32 


4C456nlln32  +  8C456n21n13  * 


(28f) 


'44n13  +  2C144nlln13  +  2C144n22n13  +  2C344n33nl  3  “  4C456n22n13 


4C456nlln13  +  SC456n21 n32  * 


(289) 
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3  (ft 

3n23 


2C44n23  +  2C144nlln23  4  2C144n22n23  +  2C344n33n23  +  4C456n22n23 
"  4C456nlln23  +  8C456n12ri31  (28h) 


30 

3n33 


12- 

C13nll  +  C13n22  +  C33n33  +  1  C333n33  +  C1 33^11 n33  4  C133n22n33 

12  12  2  2 
+  1  C113nll  4  C113n11n22  +  2  C113n22  4  C344n23  4  C344n32 


+  C344n132  +  C344n312  +  \  C366nll2  4  7  C 


366r'22  "  C366nlln22 


+  r  2  r  2 

L366n12  L366n21 


(28i) 


To  derive  expressions  for  T. the  given  by  Eqs.  (28)  and 

ij  3nij 

the  components  of  the  Jacobian  matrix  are  required.  The  components  of 
the  Jacobian  matrix  In  terms  of  the  displacement  gradients  are  (from 
Eq.  (4)): 

J11  *  1  +  ua  *  1  +  3U/3a 

J12  *  ub  *  3U/3b 
J]3  *  Uc  -  3U/3C 

J21  *  Va  '  3V/3a 

J22  =  1  +  Vb  *  1  +  3V/3b  (29) 

J23  =  Vc  *  3V/3C 
J3I  r  Wa  *  3W/3a 

J32  *  *  3W/3b 

J33  =  1  +  W£  =  1  +  3W/3C 


The  strain  parameters  can  be  expressed  in  terms  of  displacement 
gradients  if  we  substitute  (4),  (5)  and  (6)  into  (2).  Substitution 
and  simplification  give  the  following  expressions  for  n.-.- 

'  J 

”11  *  ua  ♦  !<u.2  *  \2  +  “a2> 

"12  *  "21  =  ?<ub  *  Va  *  uaub  +  VaVb  +  Vb> 

"13  '  "31  ’  >c  +  “a  *  Ve  +  Vc  +  W  (30 

"22  *  Vb  +  ?<Ub2  +  Vb2  *  “b2> 

"23  '  "32  *  !<Vc  +  “b  +  Vc  +  Vc  *  Vc> 

"33  *  “c  +  1<UC2  +  Vc2  +  “c2> 

To  evaluate  all  the  T-.'s,  we  need  the  quadratic  terms  also  in 

*  J 

p 

n^;  f.e.,  we  have  to  evaluate  ^12 *  etc*  Retaining  terms 

only  up  to  the  quadratic,  we  get  the  following  expressions  for  these 
terms. 


2 

nn  * 

^a2 

(31a) 

2 

n22  ~ 

Vb2 

(31b) 

2  _ 
n33  ' 

"c* 

(31c) 

nlln22 

■  Vb 

Old) 

nlln33 

Vc 

Ole) 

n22n33 

’  Vc 

Olf) 

riT-j2  * 

-  t(Uk2  +  v  2  +  2W.) 

Olg) 

1  =  i<Uc  +  V  +  2W 

2 2  =  K2  +  Wb2  +  2VcWb> 
nlln12  *  l(UaUb  +  UaVa) 
nlln31  =  7(UaUc  +  W 
nlln32  =  I(UaV:  +  UaV 

n22n21  =  l^VbUb  +  VaV 

n22n31  =  ?*VbUc  +  VbWa* 
n22n23  =  l*VbVc  +  VbV 
n33n21  =  I(WcUb  +  WcV 
n33n31  =  ^wcUc  +  WaWc) 
n33n32  =  l(WcVc  +  W 


(UJL  +  M,  +  V  V  +  V  W  ) 


(UJL  +  UkWK  +  V,¥  +  V  wj 


(U  V  +  U  W.  +  W  V  +  M,W.  ) 
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The  strain  parameters  are  expressed  In  terms  of  the  displacement 
gradients  by  the  two  sets  of  Eqs.  (30)  and  Eqs.  (31).  Now,  these 
expressions  are  to  be  substituted  Into  the  components  of  the  strain 
derivatives  of  the  elastic  energy  given  by  the  set  of  Eqs.  (28).  The 
resulting  expressions  and  the  components  of  the  Jacobian  matrix  given 
by  Eqs.  (29)  are  substituted  Into  the  set  of  Eqs.  (20)  to  obtain  the 
expressions  for  the  components  of  the  stress  tensor.  After  simplifi¬ 
cation,  the  components  of  the  stress  tensor  are: 


TnsUa*Cn+VC12  +  Mc*C13 

+  Ua2(l  C11  +  7  Clll  +  7  C1 66  +  7  C266 

+  Va2(i  C11  +  7  C166  ■  7  C266) 

+  Wa2(7  C11  +  7  C144  +  C456) 

+  +  ?c166  -?c266) 

+  Vb2(7  C12  +  7  Clll  "  7  C1 66  +  7  C266 

+  wb2(^  c12  +  \  c144  -  c456) 

+  Uc2(7  C13  +  C44  +  7  C1 44  +  C456^ 

+  Vc^7  C13  +  7  C144  "  C456^ 

+  Wc2(7  C13  +  7  C1 33^ 

+  UaVb^C12  +  Clll  '  C1 66  "  C266^ 


) 


) 


+  UaWc(C13  “  C1 1 3  +  C366^ 
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. 


t  r 


I 


I- 


* 


+  VbVC144  ‘  2C456^  +  WC33  +  C344^ 


+  UbVC44  +  2C456^ 

+  UbWb(C44  +  2C456>  +  VaWb  ‘  2C456 

T21  =  Ub(l  C11  '  1  C12)  +  Va  '  C11  "  7  C12J 

+  UaVa(C1T  +  C166  ‘  C266^ 

+  UaUb*?  C11  ‘  1  C12  +  C1 66  "  C266* 

+  VaVb(C11  +  C166  +  C266^ 

+  VaWc(C13  +  C366> 

+  WaWb(?  C11  '  1  C12  +  2C4S6) 

+  UbV?  C11  "  7  C12  *  C166  +  C266^ 


+  UbWc  *  C366  +  UcWb  *  2C456 


+  UcVc(C44  +  2C456> 
+  WC44  +  2C456> 


T22  =  Ua  ’  C12  +  Vb  *  C11  +  Wc  *  C13 

+  Ua2(?  C12  +  7  CH1  ’  7  C)66  “  C266] 
+  UbZ(l  C11  +  7  C166  +  7  C266^ 

+  Uc2(?  C13  +  7  C1 44  -  C456) 


(32c) 


(32d) 


+  VbWb(Cll  +  C144  +  2C456- 
+  UaWb*C12  +  C1 44  ‘  2C456) 
+  WbWc(C13  +  2C44  +  C344> 

+  UbUc(C44  +  2C456> 

♦  Vc<C44  +  C1 44  +  2C456^ 
+  Vc(C144  '  2C456^ 

+  VcWc(C44  +  C344^ 


+  VaUc  '  2C456 


ua  *  C13  +  Vb  *  C13  +  Wc  '  C33 

+  Ua2(l  C13  +  I  C113  +  7  C366) 
+  ub2(r  C13  +  l  C366> 

+  UC2{7  C33  +  1  C344) 

+  Va2{?  C13  +  \  C366) 

*  Vb2{?  C13  +  \  C113  +  7  C366^ 

+  VC2(?  C33  +  7  C344^ 

+  waz(c44  ♦  l  c13  +  j  c344) 

+  Wb2(C44  +  ?  C]3  +  l  C344) 


+  Wc  *2  C33  +  1  C333'  +  UcWa^C44  +  C344* 

+  VcWb^C44  +  C344^  +  UaWc^C13  +  C133^ 

+  VbWc^C13  +  C1 33 ^  +  UaVC113  '  C366^  +  UbVC366*  (32 

Longitudinal  Wave  Propagation  along  the  [100]  Direction  (or  a-Axis 
We  intend  to  propagate  only  longitudinal  waves  for  two  compelling 
reasons.  First,  the  capacitive  receiver  used  in  the  measurements  responds 
only  to  longitudinal  displacements;  second,  the  transverse  modes  are 
always  coupled  to  the  longitudinal  modes  as  is  evident  from  the 
expressions  for  the  stress  tensor  given  by  (32).  If  we  calculate  the 
transverse  displacement  derivative  of  the  given  by  (32)  along  any 
direction,  we  find  that  in  the  nonlinear  terms  the  longitudinal  compo¬ 
nents  always  contribute  to  transverse  waves;  i.e.,  transverse  waves  never 

propagate  as  pure  modes.  This  is  in  agreement  with  the  theory  of 
33 

Goldberg  on  isotropic  solids.  On  the  other  hand,  longitudinal  waves 
propagate  as  pure  modes  in  certain  directions  and  are  coupled  only  to 
their  second  harmonics  in  these  directions. 

For  plane  longitudinal  waves  propagating  in  the  [100]  direction, 
the  displacement  components  become 

U  =  U(a,t),  V  =  V(a,t),  W  *  W(a,t) 

or  Ik,  U  ,  V.  ,  V.  W.  and  W„  vanish.  So  for  this  case  the  set  of 
b  c  b  c  b  c 

equations  (32)  simplifies  as  follows. 


T11  =  CllUa  +  (7  C11  +  7  C111  +  7  C166  +  7  C266)Ua 
+  (]■  cu  +  \  c166  -  \  c26g)va2 

+  (]■  Cn  +  \  C144  +  C456)Wa2  (33a) 

T12  =  7(C11  ‘  C12^Va  +  7*C11  "  C12  +  2C166  "  2C266^UaVa  (33b) 


T13  =  C44  ’  Wa  +  ^C44  +  C144  +  2C456^UaWa 
T21  “  7(C11  ‘  C12)Va  +  (C11  +  C166  "  C266)UaVa 

T22  a  C12Ua  +  (7  C12  +  7  Clll  "  7  C166  "  C266)UaZ 

+  (7  C11  +  7  C1 66  +  7  C266)Va2 
+  (1  c12  +  \  c144  -  c456)wa2 

T23  =  (C44  +  2C456}VaWa 
T31  *  C44Wa  +  *C11  +  C1 44  +  2C456*UaWa 
T32s  <IC11  '7C12  +  2WVaWa 
T33BC13Ua  +  (7C13  +  7C113  +  7C366)Ua2 

+  ^7  C13  +  7  C366^Va2 
+  ^C44  +  7  C13  +  7  C344^Wa2 


(33c) 

(33d) 


(33e) 

(33f) 

(33g) 

(33h) 


(331] 


The  equations  of  motion  for  plane  wave  propagation  In  the  [100] 


direction  are 


(34) 


PQ  U  *  31^/38 
p0  V  -  3T2<|/3a 
P0  w  *  3T31/3a 


Performing  the  differentiation  and  substituting  for  3T^/3a,  etc.,  we 
get 


+  (3C11  +  Cni  +  3C166  +  C266)UaUaa 
+  *C11  +  C1 66  "  C266)VaVaa 

+  <C11  +C144  +  2WWaWaa  (35a) 

(,o  »  ■  !<cn  -  c12)vaa  *  (Cn  *  c166  -  c266)Cuavaa*vauaa) 

(35b) 

p0  5  =  C44  '  Waa  +  <C11  +  C144  +  2C456KUaWaa  +  WaUaa>  <35c> 

where 


32u 

I? 


u.  *  S 


aa 


cMJ 

3a 
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etc. 

etc. 

etc. 


For  longitudinal  waves  along  the  a-axls  we  have 
V  •  0  and  W  =  0. 

So  the  equation  of  motion  for  this  case  reduces  to 


T  'JU11  T  ''111  T  ^166  T  ^266'  aa  7T 

3a 


2  2 
3  U  3  U  .  . 


2 

3U  3  U 

3a  ^7 


“  =  C11 

6  *  3C|i  +  Cni  +  3C]66  +  C266 


B.  Longitudinal  Wave  Propagation  along  the  [010]  Direction 

For  this  case  the  only  nonvanishing  displacement  components  are 

U  *  U(b,t),  V  *  V(b,t),  W  =  W(b.t)  . 

So  for  this  case  the  T^'s  given  by  Eqs.  (32)  simplify  to  the  following. 


T11 

■  C12vb 

+  (IC11  +1 

C166  ' 

7' 

:266)Ub2 

♦  <7 

ci2 + 1  Cm  " 

7  C1 66  + 

1  C266^Vb2 

*<7 

C12  +  \  C144  " 

C456> 

w  2 

wb 

(38a) 

T12 

“  r<cii 

*  C12)ub  +  (C 

n  +  c 

166 

+  C266^UbVb 

(38b) 

T13 

■  <C44 

+  2C456>UbWb 

(38c) 

T21 

■  7(cll 

"C12^Ub+  ^7  C 

i 

n  "2 

C12 

+  C166  +  C266)UbVb 

(38d) 

T22  * 

cllvb  + 

(?  cn 

+  1  C1 66  +  1  C266)Ub 

+ 

(f  C11  + 

1  cni 

+  f  C166  ”  \  C266)Vb2 

+ 

(lcn  + 

1  C144 

+  C456)Wb2 

T23  ’ 

C44  ‘  Wb 

+  (c44 

+  C144  +  2C456)VbWb 

T31  * 

(lcn  * 

\  C12 

+  2C456^UbWb 

II 

CVJ 

CO 

h- 

C44Wb  + 

(cn  + 

C144  +  2C456^VbWb 

T33* 

C13Vb  + 

C13 

+  \  C366)Ub2 

+  C13 

+  rC113  +  ?C366)Vb2 

+  (!  C13 

+  C44 

+  C344,Wb2  * 

The  equations  of  motion  for 

plane  wave  propagation  in  the  [010] 

direction  are 

pU  * 

3T12/3b 

pii  * 

3T22/3b, 

and 

pw  * 

3T32/3b  . 

Differentiating  (38b) 

,  (38e) 

,  and  (38h)  we  get 

3T12 

3b 

c12) 

Ubb 

♦  (c„ 

+  C->c£)(UhVhh  +  VJJhh) 
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~W  s  CllVbb  +  (Cn  +  C1 66  +  C266^UbUbb 
+  ( 3C-j  i  +  Cin  +  3C166  -  C26g)VbVbb 

+  (C^  +  C144  +  2C45g)WbWbb  (40b) 


3T32  .  f 
3b  ‘  l44 


Wbb  + 


(C 


11 


+  C'-i,  +  2C456^VbWbb  +  WbVbb^ 


'144 


(40c) 


For  longitudinal  waves  In  this  direction,  we  have  11=0  and  W  =  0. 
So  the  equations  of  motion  reduce  to 


or 


r  32V  _  c 


32V 
11  3? 


+  <3C11  +  Cni  +  3C166  -  C266) 


'111 


3V 

3b 


32V 

17 


32V 

7““ 


32V 

^7 


6  Tb 


32V 

"a? 


where 

a  =  Cn 
and 


$  s  3C,i  ^  Cm  ^  3C  1  CC  "  C 


'11 


'111 


166  266 


(41) 


(42) 


C.  Longitudinal  Wave  Propagation  along  the  [001]  Direction  (or  the 
C-Axis) 

For  this  special  case  the  only  nonvanishing  displacement 
components  are 


U  *  U(c,t),  V  =  V(c,t),  W  =  W(c,t) 


J  ,  U.  ,  V,,  V. ,  W.  and  WK  vanish.  To  the  T.  .'s  for  this  case  reduce 
a  b  aba  b  ij 


to 


T11  =  C13  *  Wc  +  C13  +  C44  +  1  C144  +  C456)Uc2 


+  ^C13  +  rC144-WVc2+^C13  +  lCl33>Wc2  <4 


T12  =  (C44  +  2C456^UcVc 


T13  =  C44  ’  Uc  +  (C33  +  C344)UcWc 


T21  =  (C44  +  2C456)UcVc 


T22  *  C13  *  Wc  +  C13  +  ?  C1 44  "  C456)Uc2 

+  ci3  *  c44  +  I  C144  +  C456)Vc2 
+  4  C13  +  I  CT33^WC2 


T23  =  C44Vc  +  (C33  +  C344>VcWc 


T31  =  C44Uc  +  (C44  +  C344^UcWc 


(4: 

(4 

(4 


(4 

(4 

(4 


T32  =  C44Vc  +  (C44  +  C344)VcWc 

T33  =  C33Wc  +  ^  C33  +  !  C344^Uc2 

+  C33  +  \  C344^Vc2  +  C33  +  \ 


(4 


(4 


The  equations  of  motion  for  plane  wave  propagation  In  the  [001] 
direction  are 
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°  ’  C33 
and 

6  ~  3C33  +  C333  ' 


D.  Solution  to  the  Equations  of  Motion 

The  equations  of  motion  for  longitudinal  waves  propagating  along 
the  [100],  [010],  and  [001]  directions  given  by  Eqs.  (37),  (42),  and 
(48)  have  the  same  general  form  given  by 


Pn  U  -  a  U  -  «  Oa 


So  we  see  that  a  pure  mode  longitudinal  wave  may  propagate  in  a  non¬ 
linear  hexagonal  solid  along  the  three  directions  [100],  [010],  and 
[001].  As  other  workers  have  noticed  for  cubic  crystals  [20],  the  non¬ 
linear  term  causes  the  wave  to  generate  harmonics  as  it  propagates.  In 
order  to  solve  Eq.  (49),  a  perturbation  approach  Is  necessary.  Let  us 
try  the  solution 

U  =  U]  +  U2  .  (50) 

Here  U-j  and  Ug  represent,  respectively,  the  approximate  solutions  of 
first  and  second  order  of  the  equation.  By  substituting  Eq.  (50)  into 
Eq.  (49)  one  gets 


"0  U1 


°0  U2  ■  °  Ulaa 


-  °  u2aa  *  5(ula 


+  U2a>(Ulaa 


*  u2aa>  <51> 


Keeping  only  first  order  terms  leads  to  the  linear  wave  equation: 
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p0  “l  *  «  Ulaa  =  0  *  <52> 

where  (a/p^)^2  *  Cq  is  the  phase  velocity  of  the  linear  wave.  Sub¬ 
stituting  (52)  into  (51),  we  have 

p0  V  “  U2aa  =  ^UlaUlaa  +  (U2aUlaa+UlaU2aa)+U2aU2aa3  <53> 

Since  U2  «  Uj,  the  second  and  third  terms  on  the  right  hand  side  are 
very  small  compared  to  the  first  term.  Thus,  to  a  second  approximation 
the  solution  is  reduced  to 


p0  U2  -  “  U2aa  *  6  UlaUlaa  * 


We  Impose  the  boundary  condition  that 


Uj  *  A  sin  wt 

u2  -  0 


where  a  =  0 


(54) 


(55) 


Thus,  a  solution  for  the  linear  equation  (52)  which  satisfies  the 
boundary  condition  Is 

U1  =  A  sin(ka  -  at)  .  (56) 

Now,  substituting  (56)  Into  (54),  we  may  rewrite  as  follows: 

p0  U2  -  a  U2ag  =  -  \  6  k3A2  sin  2(ka  -  tut)  .  (57) 

Taking  into  account  the  boundary  condition  (55)  for  U2,  we  use  a  trial 
solution  in  the  following  form: 
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i  a 

i  r 

»’ 

r 


i  r 

!  ■ 

* 

»*  - 

i* 

I  ■ 


*  i- 


U2  =  B  a  sin  2(ka  -  ut)  +  C  a  cos  2(ka  -  wt)  (58) 

where  B  and  C  are  coefficients  to  be  determined.  Putting  the  solution 
into  (57),  we  get 

4(-u^pQBa  +  k2aBa  +  kaC)sin  2(ka  -  wt)  + 

3  2 

4(-u>2p0Ca  -  kaB  +  k2aCA)cos  2(ka  -  ut)  ■  •  sir,2(ka  -  <*>t )  . 


Now,  equating  the  coefficients  on  both  sides  of  the  equation  we  get 


-  u)2p0Ba  +  k2aBa  +  kaC  =  -  ~  6k3A2' 
and 

2  2 
-oj  pqC  -  kaB  +  k  aCa  =  0 

But 


and 

to  =  Cgk  , 

so  Eqs.  (59)  become 

-  C02k2pQBa  +  k2pQC02Ba  +  kP()C02C  =  -  |  k3A2 
and 

CQ2k2p0Ca  -  kP()C02B  +  k2p0C02Ca  =  0  . 
Dividing  Eq.  (59a)  by  kpgCg2,  we  get 

c-iMi 

fip0C0 


(59) 


(59a) 


(59b) 


(60a) 
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A1 


so,  dividing  Eq.  (59b)  by  kpCQ  ,  we  get 


B  *  0  . 


(60b) 


So,  after  one  iteration  the  solution  to  the  equation  of  motion 
becomes 

U(a,t)  *  A  sin(ka  -  <*>t) 


2 

-  cos  2{ka  -  «t) 


6p0C0 


(61) 


The  equation  of  wave  propagation  in  the  nonlinear  regime  as 
given  in  Eq.  (49)  was 


d2u  d^u  _  .  du  d2u 
,20 


(62) 


Breazeale  and  Ford  write  the  same  equation  from  the  fluid  analogy 
(for  cubic  crystals)  as 


d  u  _  „  /d£u  ,  ,  du  deu\  .  v  du  d^u 

‘»^'K2<57  3^'^’  k3^  3? 


/d2u 


d2u, 


d2u 


or 


2  2  2 
d  u  y  d  u  _  Mu  «  y  \  du  d  u 

(3K2  *  K3>  3? 


(63) 


Equations  (62)  and  (63)  are  the  same  if  we  put 


a  «  K, 


and 


6  *  (3K2  +  K3) 


(64) 


Now,  If  we  consider  the  linear  wave  equation  (Eq.  52)  we  can  write  it 


E.  The  Ultrasonic  Nonlinearity  Parameter  for 
Hexagonal  Crystals 

As  in  the  case  of  cubic  crystals  [1],  we  can  define  the 
ultrasonic  nonlinearity  parameter  for  hexagonal  crystals  as  the  nega 
tive  of  the  ratio  of  the  nonlinear  term  to  the  linear  term  in  the 
nonlinear  wave  equation  [Eq.  (63)3.  It  is  given  by 


6  *  -  £  .  (69) 

For  an  initially  sinusoidal  disturbance  at  a  *  0,  the  solutions  as 
given  by  Eq.  (68)  can  be  written  in  terms  of  the  nonlinearity  parameter 
as 


A,2k2a 


1  ° 

u  «  A-j  s1n(ka  -  <nt)  +  — g —  B  cos  2(ka  -  o>t) 


(70) 


A,  k  a 

where  A-j  Is  the  amplitude  of  the  fundamental  wave  and  A?  *  — —  B 
is  the  amplitude  of  the  generated  second  harmonic  wave.  In  terms  of  A-j 
and  ?2»  *s  9iven  by 


-  «<4>  -f 

k*a 


(71) 


So  by  measuring  A^  and  Ag  one  can  determine  B  which  can  be  used  to 
evaluate  the  6's  which  are  combinations  of  SOE  and  TOE  constants,  or 
the  K3's  which  are  combinations  of  TOE  constants  only.  The  parameters 
a  (or  Kg),  6  (or  (3Kg  +  K^))  and  K^  (or  (6  -  3Kg))  for  the  three 
symmetry  directions  considered  in  Sections  A,  B,  and  C  given  by  Eqs. 
(37),  (42),  and  (48)  are  written  In  a  tabular  form  In  Table  I.  Note 


that  and  Cggg  are  combinations  of  TOE  constants  (given  by 


The  K2  and  K3  Parameters  for  Hexagonal  Crystals  along  the  Directions  Considered  In 
Sections  A,  8,  and  C 


Eqs.  (26)).  In  Table  I  we  have  given  the  Kj  parameters  in  terms  of 
these  combinations  as  well  as  In  terms  of  the  independent  TOE 
constants. 

F.  The  Equations  of  Motion  when  the  Coordinate  System 
is  Rotated  by  60°  about  the  c-axis 

In  the  remaining  calculations  we  will  consider  the  nonlinear 
equation  of  motion  for  various  directions  in  the  basal  plane.  As 
was  shown  by  Borgnis6  and  Brugger,7  any  direction  in  the  basal  plane 
is  a  pure  mode  direction  In  the  linear  approximation.  By  determining 
the  nonlinear  equation  of  motion  for  several  directions  we  can  gain 
insight  about  the  way  the  nonlinearity  parameter  changes  as  the 
propagation  direction  changes  In  the  basal  plane. 

In  order  to  derive  the  equations  of  motion  for  waves  generating 
along  the  a'-axis,  we  have  to  derive  the  expression  for  strain  energy 
in  the  frame  of  reference  in  which  the  a-axls  in  the  original  frame 
Is  rotated  by  60°  about  the  c-axis.  Consider  the  propagation  of 
plane  finite  amplitude  waves  along  the  a'-axis.  Rotate  the  basal  plane 
In  the  original  frame  by  60°  or  about  the  c-axis.  So  the  coordinate 
frame  formation  can  be  written  as 


a' 

a 

b' 

*  (R) 

b 

c' 

c 

where  the  rotation  matrix  Is  given  by 


where  8  is  the  angle  through  which  the  original  basal  plane  has  been 
rotated  about  the  c-axis.  We  rotate  the  a-axis  by  e  =  60°  in  which 
case 


(R) 


1/2  /37?  0 

-  yj/i  i/2  o 

0  0  1 


(74 


Let  us  define  the  displacement  components  in  the  new  coordinate  system 
as 


u'  *  u'  (a' ,t) 

v'  =  v'  (a',t)  (75 

w'  *  w'  (a',t) 

etc. 

Using  the  transformation  matrix  as  given  by  Eq.  (74)  one  can 
transform  the  strain  components  by  [9] 


(n)  =  (R*)(n')(R)  , 


(76 


to  arrive  at  the  following  expression: 


1/2 

/S/2 

0 


-/S/2 

1/2 

0 


1/2 

-/S/2 

0 


/3/2  0 

1/2  0 
0  1 


(78 


Multiplication  and  simplification  lead  to  the  following  expressions 
for  the  strain  components  in  the  original  frame  of  reference  expressed 
in  terms  of  the  strain  components  in  the  rotated  or  primed  coordinate 
system. 


1  i  J5  i  /S  •  .3  • 

nll  ~  T  nll  4  n12  “  7  n21  7  n22 

_  •  .  1  •  3  «  * 

*12  "  ?  *11  T  *12  ~  7  n21  ‘  7  *22 

1  •  /S  * 

n13  2  *13  "  7  *23 

_  /S  »  3  •  .  1  »  i 

n21  "  7  *11  '  4  *12  7  *21  7  n22 

_3i,/S',/S«l  • 

n22  "  7  *11  7  *12  7  *2 1  7  *22 

_  /s  •  .  1  • 

n23  ?  ”13  r  *23 

_  1  i  /S' 

*31  "  I  *31  ‘  I  *32 

_  /S  •  .  1  i 

*32  '  7  *31  7  *32 


(79 


Substitution  of  Eq.  (79)  in  Eq.  (25)  gives  the  following  expression 
for  the  strain  energy  in  the  primed  coordinate  system. 


4»*  =  ?  Cn{^ll )2  +  ^22^  +  ^n12^2  +  ^n21^^ 


+  CioCn^n 


7  (n(9)2  -  \  (nAi)2} 


12|T1lln22  ‘  I  'n12'  "  2  'n21 


+  C13{nlln33  +  T)22T)33}  +  \  C33^n33^ 


+  C44  {^n13^  +  ^31^  +  ^n23^  +  ^32^} 

+  F  Clll^nll^3  +  3(nlT^2r>22  +  3(rl22^2r,ll  +  ^22^^ 

+  7  C113n33{(nll)2  +  2nlln22  +  (n22)2} 

+  7  Cl33^n33^  ^11  +  ^33^  n22^ 

+  C144^t,11  +  +  ^n32^2  +  ^n13^2  +  ^n31  ^  } 

+  7  C166{^11^3  +  ^n22^3  +  2(n^2)2n^1 

+  2(n^)2n^  -  (n^)2^!  -  (^i)2^ 

+  2(n]2)2Ti22  +  2(n^  J2^} 

»lp  r(  i  »3  ,  3<^5‘  /  i  \3  ,  /  i  \3 

+  5C266{(nll)  +  '7_(ti12)  +  ~F(ti21) 

i  i  \3  3  /  i  \2  i  3  /  i  \2  i 

■  vn22 7  “  j  *T112'  ^11  "  7  '^21  ^11 


&  ^ 

1  i  m 

— 

1 

: 

48 

\  Q 

■*■  3 ( t^22 )  nll  ”  3(n^)  ^2  ”  9rinT,12T121 

-*  —  — 

> 

3/5  ,  i  »2  i  3/5  /  i  \2  , 

~  2  '^21 1  n12  ”  2  '^12'  ^21 

1 , 

*  2"  (n^-j )  ^22  *  5-(n{2)  r'22  2^21  ^22 ^ 

— - 

i 

I  - 

+  S  C333^n33^ 

•  *  « 

1 

+  C344T133^T123^  +  ^32^  +  ^13^  +  ^n31^  * 

\  . 

i  r 

+  1  C366n33^nll^  +  2^i2^  +  ZU^-j  )2  +  (n22)  " 

2nlln22} 

L  i 

» 

r* 

+  \  C456{[1  (nil }  *  T  n12  +  T  n21  *  ?  n22] 

;  • 

IK 

x  f?  ^i3)2  "  \  (n23)2  +  ^  T'73n23  +  ?  (n37'2 

l 

-  ^  (n^)  +  ^ 
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Once  the  expression  for  strain  energy  is  obtained,  one 

can  execute 

the  algebra  exactly  as  in  Sections  A,  B,  and  C  and  derive  the  equations 
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of  motion  for  plane  pure  mode  longitudinal  waves  propagating  along  the 
[100],  [010],  and  [001]  symmetry  axes  in  the  rotated  frame  of  reference. 
This  has  been  done  and  we  found  that  for  these  directions  also  the 
equation  of  motion  has  the  same  form  as  Eq.  (49)  and  its  solution 
follows  exactly  as  in  Section  D.  The  parameters  a  and  6  appearing  in 
the  equation  of  motion  and  the  corresponding  Kg  and  parameters  are 
tabulated  in  Table  II. 

From  Tables  I  and  II  one  can  notice  that  the  a-axis  and  a1 -axis 
or  the  [100]  direction  in  the  original  and  rotated  coordinate  systems 
are  equivalent.  Similarly,  the  directions  perpendicular  to  them  which 
are  the  [010]  directions  in  the  two  cases  are  also  equivalent.  Of 
course,  the  c-axis  is  the  same  as  before;  e.g.,  the  rotation  operation 
is  performed  about  the  c-axis. 


6.  The  Equations  of  Motion  when  the  Coordinate  System 
is  Rotated  by  -it/4  about  the  c-Axis 

In  order  to  derive  the  equations  of  motion  for  the  propagation  of 
plane  longitudinal  waves  along  the  [110]  direction  in  the  basal  plane, 
we  have  to  derive  the  expression  for  strain  energy  in  the  frame  of 
reference  in  which  the  a-axis  in  the  original  coordinate  system  is 
rotated  by  45°  about  the  c-axis.  The  coordinate  transformation  for  this 
case  is: 


a1 

a 

b' 

=  (R) 

b 

c' 

c 

I 


Table  II.  The  Kz  and  K3  Parameters  for  Hexagonal  Crystals  In  the  Rotated  Frame  of  Reference 
Described  in  Section  F 


where 
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MfZ  Vfl  0 

R  *  -7 !/[  ]/S?  0  (82) 

0  0  7 

is  the  transformation  matrix  from  Eq.  (73).  The  components  of  strain 
in  the  initiai  frame  can  be  expressed  in  terms  of  the  components  of 
strain  in  the  rotated  frame  as 

n  -  R*  (n‘)  R  (83) 


n17  n72  n73|  -7/*^  0  7//Z  1/^2  0 

n27  n27  n23  =  "\/^2  0  n 22  n23  " M^2  0 


n37  n32  n33 


0  1  n37  n32  n33  0  00 


Expansion  and  simplification  gives  the  foilowing  expressions  for  n^j's 


in  terms  of  n^'s. 

*  J 
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*  n23^ 

n21 

=  \  (n^  - 

n72  + 

n27 

’  n22* 

n22 

=  \  ^71  + 

n72  + 

n27 

+  n22) 
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n23  =  ^  (n13  +  n23^ 
n31  =  ^  *n31  *  ^32  ^ 
n32  *  ^  (r,31  +  n32J 


n33  *  n33 


Substituting  Eq.  (85)  into  Eq.  (25),  we  get  the  following  expression 
for  the  elastic  strain  energy  in  the  rotated  or  primed  frame  of 
reference. 


-  1  r  r_i  2.1  2.i  2.1  2-, 

^  CllLnll  n22  n12  n21  J 


+  r  ( n 1  n 1  1  I  ^  1  1  2-1 

C12$nlln22  ?  n12  *  ?  n21  ] 


+  C13^nlln33  +  n22n33^  +  I  C33n332 
+  C44[n132  +  n312  +  n232  +  n322]  +  1  C333n333 
+  1  Clll^ll3  +  3nll2  n22  +  3rt222nll  +  n223^ 


+  2  C133n332^nll  +  n22^ 


+  I  C113  n33^il2  +  2niln22  +  n22^ 

+  C144^nl  1  +  n22^n232  +  n322  +  n132  +  n31^ 
+  C344n33^n232  +  n322  +  n132  +  n312^ 
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The  derivatives  of  the  strain  energy  with  respect  to  the  strain 
components  are  obtained  by  taking  the  approximate  derivatives  of 
Eq.  (86).  Then,  substituting  them  in  Eqs.  (20),  one  obtains  the  com¬ 
ponents  of  stress  tensor.  These  stress  components  are  differentiated 
with  respect  to  the  appropriate  displacement  to  obtain  the  equations 
of  motion  for  longitudinal  waves  propagating  along  [100],  [010],  and 
[001]  directions  in  the  rotated  frame  which  are,  respectively,  the 
[110],  [TlO],  and  [001]  directions  in  the  original  coordinate  system. 
The  algebra  is  similar  to  that  performed  in  Sections  A,  B,  and  C,  and 
is  not  reproduced  here. 

It  is  found  that  for  these  directions  also,  the  equations  of 
motion  have  the  same  form  as  Eq.  (49)  and  the  solution  follows  as  in 
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Section  D.  So  It  Is  found  that  a  pure  mode  longitudinal  waves  propagate 
along  these  directions  with  the  generation  of  harmonics.  The  parameters 
o  and  6  appearing  in  the  equations  of  motion  and  the  corresponding  Kg 
and  Kg  parameters  are  tabulated  in  Table  III. 

An  examination  of  Table  III  reveals  that  for  the  [100]  and  [010] 
directions  in  the  45°  rotated  frame,  the  Kg  and  Kg  values  are  the  same. 
This  can  be  expected  from  the  symmetry  of  hexagonal  crystals.  Obviously, 
the  values  for  the  c-axis  do  not  change,  because  the  rotation  was  about 
the  c-axis.  The  Kg  values  are  the  same,  i.e. ,  C^,  for  every  direction 
in  the  basal  plane.  This  is  consistent  with  the  fact  that  the  velocity 
of  longitudinal  waves  is  the  same  {V  *  ^C-^/pq)  for  every  direction  in 
the  basal  plane. 

H.  Summary  of  Results 

We  can  arrive  at  some  very  important  conclusions  by  comparing 
Tables  I,  II,  and  III.  They  predict  the  relative  magnitudes  of  the 
second  harmonic  generated  when  waves  are  propagated  along  different 
directions  in  a  hexagonal  crystal.  The  Kg  parameters  for  a-axis  and 
a'-axis  are  the  same,  i.e.,  (2Cggg  -  C^g).  The  Kg  values  for 
directions  perpendicular  to  them  (the  [010]  direction)  are  also  the 
same,  i.e.,  (5/2  Cggg  -  C^-j  -  1/2  C^g).  These  directions  make  a  30° 
angle  with  the  a-  and  b-axes.  By  symmetry  the  30°,  90°,  and  150° 
directions  are  equivalent  and  they  have  the  same  Kg  parameter.  The  Kg 
values  for  the  [110]  direction  and  [T10]  directions  which  make  a  45° 
angle  with  the  a-axis  (and  15°  with  the  b-axis),  as  given  in  Table  III, 


Described  In  Section  G 


are  the  same,  i.e.,  ^  (9C222  "  3C112  ”  2Cill^  By  symmetry  the 
directions  which  make  15°,  45°,  75°,  105°,  etc.  are  equivalent  and 
they  have  the  same  ICj  parameter.  Taking  into  account  the  symmetry 
of  different  directions  we  have  combined  Tables  1,  II,  and  III  to 
make  Table  IV,  which  is  a  summary  of  the  results. 

Another  Important  observation  that  can  be  made  is  that  the 
difference  in  K3  parameters  between  that  for  the  a -axis  and  that  for 
the  45°  direction  is  the  same  as  the  difference  between  the  K3 
parameters  for  the  45°  direction  (or  15°  direction)  and  the  30° 
direction.  The  value  of  the  difference  is  C266  or  ^  (2C^ ^ ^  -  C112 
-  C222).  Note  that  always  the  angles  of  the  various  directions  are 

measured  from  the  a-axls  which  is  taken  as  the  [100]  direction.  The 
average  of  the  K3  parameters  for  [100]  and  30°  direction  is  the  K3 
parameter  for  15°  direction  (or  45°  direction).  This  suggests 

5 

a  diagram  similar  to  the  slowness  diagrams  of  Musgrave,  as  shown  in 
Figure  2.  However,  in  Figure  2  we  Indicate  the  magnitude  of  the  K3 
parameter  in  the  basal  plane.  The  slowness  diagram  for  longitudinal 
waves  in  the  basal  plane  of  a  hexagonal  crystal  would  be  a  circle 
because  the  velocity  is  the  same  In  all  directions.  Even  though  the 
velocity  is  the  same  In  all  directions,  the  amplitude  of  the  generated 
second  harmonic  varies  with  direction,  as  indicated  in  Figure  2.  This 
means  that  by  measuring  the  second  harmonic  and  calculating  K3  along 
the  [100],  the  [110]  and  the  [010]  directions  in  a  hexagonal  crystal 
one  should  be  able  to  Isolate  C266  and  the  combination  (C-j^  +  3C366) • 
or  equivalently,  the  combinations  {^2ZZ  ~  C112^  and  \  ^3C112  "  *C111^' 
Measuring  K3  along  the  c-axis  should  give  C333- 


Table  IV.  The  K2  and  K3  Parameters  for  Hexagonal  Crystals  (Summary  of  Tables  I,  II,  and  III) 
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Finally,  according  to  our  Interpretation  of  the  figure  giving 

the  pure  mode  directions  In  a  hexagonal  crystal  In  the  publication  of 
15 

Brugger,  there  remains  a  set  of  pure  mode  directions  for  which  we 
have  not  yet  been  able  to  derive  a  nonlinear  theory.  These  are  the 
directions  along  a  cone  whose  apex  angle  Is  centered  about  the  c-axis. 
Such  a  derivation  may  be  contained  In  a  later  technical  report;  how¬ 
ever,  at  present  It  appears  that  the  theory  may  become  so  complicated 
that  one  can  derive  the  results  only  for  specific  crystals  rather 
than  In  general  terms  as  given  above  for  the  directions  we  have 
considered. 
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